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Abstract — The maximal rate for non-square Complex Orthog- 
onal Designs (CODs) with n transmit antennas is i + ^ if n is 
even and \ + ^-j- if n is odd, which are close to \ for large 
values of n. A class of maximal rate non-square CODs have been 
constructed by Liang (IEEE Trans. Inform. Theory, 2003) and 
Lu et. al. (ffiEE Trans. Inform. Theory, 2005) have shown that 
the decoding delay of the codes given by Liang, can be reduced 
by 50% when number of transmit antennas is a multiple of 4. 
Adams et. al. (IEEE Trans. Inform. Theory, 2007) have shown 
that the designs of Liang are of minimal-delay for n equal to 
1 and 3 modulo 4 and that of Lu et.al. are of minimal delay 
when n is a multiple of 4. However, these minimal delays are 
large compared to the delays of the rate 1/2 non-square CODs 
constructed by Tarokh et al (IEEE Trans. Inform. Theory, 1999) 
from rate-1 real orthogonal designs (RODs). In this paper, we 
construct a class of rate-1/2 non-square CODs for any n with 
the decoding delay equal to 50% of that of the delay of the 
rate-1/2 codes given by Tarokh et al. This is achieved by giving 
first a general construction of rate-1 square Real Orthogonal 
Designs (RODs) which includes as special cases the well known 
constructions of Adams, Lax and Phillips and Geramita and 
Pullman, and then making use of it to obtain the desired rate- 1 
non-square COD. For the case of 9 transmit antennas, our rate-| 
COD is shown to be of minimal-delay. The proposed construction 
results in designs with zero entries which may have high Peak-to- 
Average Power Ratio (PAPR) and it is shown that by appropriate 
postmultiplication, a design with no zero entries can be obtained 
with no change in the code parameters. 



I. Introduction and Preliminaries 

There are several definitions of Orthogonal Designs (ODs) 
in the literature (TJ, Q, O the well known being as given in 

ID: 

Definition 1: A Complex Orthogonal Design (COD) 
G(xo,X\, ■ ■ ■ ,Xk-i) (in short G) for n transmit antennas is 
defined as a p x n matrix such that (i) the nonzero entries of 
G are the complex variables ±xo,±xi, ...,±Xk-i and their 
conjugates and (ii) G n G — (|xo| + • ■ • + |a:fc_i| )In where 
H stands for the complex conjugate transpose and I n is the 
nxn identity matrix. The matrix G is also said to be a [p, n, k] 
COD and its rate in complex symbols per channel use is -. 
When xo, ■ ■ ■ ,Xk-i are real variables, the designs are called 
Real Orthogonal Design (ROD). 
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Space-time block codes (STBCs) from CODs have been 
widely studied for square designs, i.e., p = n, since they 
correspond to minimum decoding delay codes for co-located 
multiple-antenna coherent communication systems. However, 
non-square designs naturally appear and important in the 
following situations. 

1) In coherent co-located MIMO systems, for a specified 
number of transmit antennas, non-square designs can 
give much higher rate than the square designs 0] . 

2) In non-coherent MIMO systems with non-differential 
detection, non-square designs with p = 2n lead to low 
decoding complexity STBCs J2). 

3) Space-Time-Frequency codes can be viewed as non- 
square designs J3]. 

4) In distributed space-time coding for relay channels, 
rectangular designs naturally appear |4|. 

5) Rate | non-square CODs have been proposed for use 
in analog transmission with application to channel feed- 
back 0. 

The rate of the square CODs falls exponentially with 
increase in the number of transmit antennas. The following 
theorem relates the rate of a square OD, real/complex, with the 
number of transmit antennas. For an integer n = 2 a (2b + 1), 
where a = 4c + d and a,b,c and d being integers with 
< d < 3, the Hurwitz-Radon number p(n) is defined as 
p(n) = 8c + 2 d . 

Theorem 1 ( f@j, fiT]/, ^j): The maximal rate of a square 
ROD for n transmit antennas is given by ^^- where p(n) is 
the Hurwitz-Radon number of n, while that of a square COD 
is given by ^ii where a is the exponent of 2 in the prime 
factorization of n. 

Several authors have constructed square CODs achieving 
maximal rate [6|, [§). In (6), the following induction method is 
used to construct square CODs for 2 a antennas, a = 2, 3, • • • , 
starting from 
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where G a is a 2 a x 2 a complex matrix. Note that G a is a 
square COD in (a + 1) complex variables xq, X\, X2, • ■ • , x a . 
It is clear from Theorem Q] as well as the construction 
given by (Q]) that the square ODs, real and complex are not 
bandwidth efficient and naturally one is led to study non- 
square ODs in order to obtain codes with higher rate. It is 
known that [7| there always exists a rate-1 ROD for any 
number of transmit antennas. In fact, the existence of rate- 
1 [p, n,p] ROD is equivalent to that of a \p,p, n] square ROD. 



For a rate-1 ROD, the minimum value of decoding delay p as 
a function of n, denoted by v(n), is given by 



v(n) 
S(n) : 
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4s if n ■■ 
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As + 3 if n : 



-1 
•2 

- 3 or 8s + 4 
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It is known [1J that the maximal rate of a non-square COD 
is equal to -| + ■£- when the number of transmit antennas is 
2a — 1 or 2a. Liang in [ 1 1 has given an explicit construction of 
non-square CODs achieving this maximal rate for any number 
of antennas. There is also another construction of these codes 
given by Lu et al iflOl . The former construction is algorithmic 
in nature while the latter one is based on patch-up of several 
matrices. The minimum decoding delay for the maximal rate 
non-square CODs is, in general, not known. The following 
theorem states what is known about the minimum delay of 
these code. For details, see ifTTIl . 

Theorem 2 ( 4771/ ): A tight lower bound on the decoding 
delay of maximum rate non-square CODs for n antennas is 
( _i) f° r n = 2a—lorn= 2a. Moreover, if n is congruent 
to 0, 1 or 3 modulo 4, then this lower bound is achievable. If 
n is congruent to 2 modulo 4, the minimum decoding delay 
is upper bounded by 2 („_*,) ■ 

As the rate of these maximal rate codes is close to 1/2 for large 
number of antennas, it is sufficient to focus on rate 1/2 codes 
when large number of antennas is under consideration. Then, a 
natural problem to study is construction of rate 1/2 non-square 
CODs with the decoding delay as small as possible. Tarokh et 
al have given a class of rate 1/2 codes obtained from rate- 
1 RODs, which has lower delays than those of maximal-rate 
codes constructed in [ 1 1, 1 10 1 for number of transmit antennas 
more than 5. For example, for four transmit antennas, the rate 
1/2 code is 



G = 



V2 



x 


-Xi 


-x% 


-X3 


Xl 


•t'o 


x$ 


-x 2 


Xl 


-xa 


x 


Xl 


X3 


■r-2 


-Xl 


Xo 


Xq 


* 

-x\ 


—x 2 


* 



(3) 



Notice that, contrary to the definition in [7| of a COD, in 
this code, the variables appear in a column more than once and 
each entry of all the columns of the design matrix is scaled 
by -y= in order to satisfy the condition G n G = (\xo\ + ■ ■ ■ + 

1^3 1 )^4 of a COD. We call such designs scaled COD which 
is not a COD in the conventional sense as in [ 1 ] (Definition 
Q]). We define the class of scaled CODs as follows: 

Definition 2: A X-scaled complex orthogonal design, for a 
positive integer A, (A-scaled COD) G is a p x n matrix in 
k complex variables xq,X\,--- ,Xk-i such that a non-zero 
entry of the matrix is a variable or its complex conjugate, 
or the negative of these and all the entries of any subset of 
columns of the matrix is scaled by -k= satisfying the condition: 



G H G = (\xq\ + ■ ■ ■ + \xk-i\ )I n - The matrix G is also said 
to be a [p, n, k] A— scaled COD. 

Notice that a A-scaled COD with with no column scaled by 
-4= is a COD. In columns with scaling by -4= all the variables 
appear exactly A times. In this paper we consider only the case 
A = 2 and call these codes simply scaled-COD. 

Contributions of this paper: The contributions of this paper 
may be summarized as follows: 

• For the rate 1/2 scaled CODs of Q, all the columns 
are scaled by the factor -4=, which led to the reduced 
delay compared to the codes of Liang and a main result 
of this paper is that by having only a subset of the 
columns scaled by -h= further reduction in delay by 50% 
is possible. We use following notations to refer to the 
rate 1/2 CODs given by Tarokh et al [7], the maximal 
rate codes given in Lu et al ifTOi l and the codes of this 
paper. 

- L n is the maximal rate COD for n transmit antennas 
with the decoding delay as specified in the Theo- 
rem 12 

- TJC n is the rate 1/2 scaled CODs for n transmit 
antennas constructed by Tarokh et al [7|. 

- (DR) n is the rate 1/2 scaled CODs for n transmit 
antennas constructed in this paper. 

Note that as n increases, the maximal rate of L n ap- 
proaches 1/2, thus two codes L n and TJC n can be 
compared for large value of n, based on their delays. 
It is not difficult to see that the decoding delay of TJC n 
is less than that of L n for large n. We provide an explicit 
construction of rate 1/2 scaled CODs for any number of 
transmit antennas such that the decoding delay of these 
codes is v(n) when n is the number of transmit antennas, 
whereas the delay for the codes TJC n is 2u(n). The 
Table U at the top of the next page shows that for large 
values of n, but for a marginal decrease in the rate with 
respect to L n , the codes of this paper are the best codes. 

• For the case of 9 transmit antennas our rate-i code is 
shown to be of minimal-delay. 

• As a byproduct of the above mentioned construction, a 
general construction of square Real Orthogonal Designs 
(RODs) is presented which includes as special cases well 
known constructions of Adams, Lax and Phillips [9] and 
Geramita and Pullman lfl2ll . 

• Even though scaling only a subset of columns allowed 
us to decrease the delay, it is shown that such a scaling 
limits the rate of the design strictly to \ . In other words, 
the maximal rate of the scaled-CODs is | when scaling 
is present in atleast one column. 

• Zero entries in a design increase the Peak-to-Average 
Power Ratio (PAPR) in the transmitted signal and it is 
preferred not to have any zero entries in the design. 
This problem has been addressed for square designs 
in |fl~3), 03). All the known maximal rate non-square 
designs have zero entries. Our initial construction of rate- 
| scaled CODs have zero entries in the design matrix 
which will lead to higher Peak-to-Average Power Ratio 
(PAPR) in contrast to the designs TJC n . However, we 
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The comparison of maximum rate achieving codes and rate 1/2 codes 
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show that by post-multiplication of appropriate matrices, 
our construction leads to designs with no zero entries 
without change in the parameters of the design. 

The remaining part of the paper is organized as follows: In 
Section [H] we present the main result of the paper given by 
Theorem[5] Before this, construction of a new set of maximal- 
rate square RODs is given in Subsection IH-AI In Subsection 
III-BI construction of two new sets of rate-1 RODs from the 
maximal-rate square RODs of Subsection IH-AI is presented 
and in Subsection IH-CI construction of the low-delay rate- 1/2 
scaled-CODs is achieved using rate-1 RODs of the previous 
subsection. In Subsection IH-DI it is shown that the maximal 
rate for scaled-CODs is | . For the special case of 9 transmit 
antennas, in Section [TTH it is shown that our construction is 
of minimal delay. In Section IIVI we show that the codes 
discussed so far can be made to have no zero entries in 
by appropriate preprocessing without affecting the parameters 
of the design. Concluding remarks constitute Section [V] In 
Appendix iBl it is shown that the well known constructions of 
square RODs by Adams-Lax-Phillips and Geramita-Pullman 
are special cases of our construction. 

II. A Construction of Rate- 1/2 Scaled Complex 
Orthogonal Designs 

In this section, we construct a rate- 1/2 scaled CODs for any 
number of transmit antennas with 50% reduction in decoding 
delay compared to the rate- 1/2 codes constructed by Tarokh et 
al [7|. The construction of these codes is done in the following 
three steps: 

STEP 1 : Construction of a new set of square RODs (Subsec- 
tion UtaJi. 

STEP 2: Construction of two new sets of rate-1 RODs from 
the square RODs of STEP 1 (Subsection llTBl . 
STEP 3: Construction of low-delay rate-1/2 scaled CODs 
using rate-1 RODs (Subsection IH-CI ). 

Before explaining these steps, we first build up some 
preliminary results needed to describe these steps. 

Let F2 = {0, 1} be the finite field with two elements 
with addition and multiplication denoted by 61 © 62 and 6162 
for 61,62 G F2. We consider logical operations also on the 
elements of this field: 61 + 62 and 61 represent respectively 
the logical disjunction (OR) of b\ and 62 and complement or 
negation of bi, i.e., 



For any finite subset B of the set of natural numbers N, 

let a G N be the least integer such that 6 < 2 a for any 6 <G B. 
Often we identify each element of B with an element of ¥ 2 us- 
ing the following correspondence: 6 G B <-* (b a -i : ■ • • , &o) G 
F2 such that 6 = X^=o bj2*, bj G F2. The all zero vector and 
all one vector in F§ are denoted by and 1 respectively. For 
x € B, x, x c and |x| represent respectively the 2's complement 
of a: in F|, l's complement of x in ¥?, and Hamming weight 
of x. In other words, x — 2 a — x and x c = 2 a — 1 — x. 
Let x = (x a -i,-" i^o) and y = (y a -i,--' ,Vo)- Then, 
x © y, x ■ y denote the component-wise modulo-2 addition and 
component-wise multiplication (AND operation) of x and y 
respectively i.e., x © y = (a; a -i © 2/a-i, • • ■ ,%o © Vo) and 
x-y = (x a -iy a -i,--- ,x y ). 

Let Z 2 a = {0, l,--- ,2 a - 1}. For a set K C Z 2 a, we 
define K = {x \ x G K}, K c = {x c | x G K}, m®K = 
{m © a I a G K} for any m G Z 2 ° and \K\ denotes the 
number of elements in the set K. 

For any two sets A, B with B C A, the set A \ B, consists 
of those elements of A, which are not in B. For two integers 
i,j, we use the notation i = j, to indicate that i — j = 
mod 2. 

For any matrix of size n,\ xn 2 , the rows and columns of the 
matrix are labeled by the elements of {0, 1, • • • , n\ — 1} and 
{0, 1, • • • , n 2 — 1} respectively. If M is a p x n matrix in k real 
variables xo,xi,x 2 ,- ■ ■ , Xk-i, such that each non-zero entry 
of the matrix is Xi or — Xi for some i G {0, 1, • • • , k — 1}, 
obviously, it is not necessary that M is a ROD. For example, 

is not a ROD. In the following (Lemma [TJi, we 
Xi x \ 

derive a necessary and sufficient condition for the matrix M 
to be a ROD. 

A submatrix M 2 of size 2x2, constructed by choosing any 
two rows and any two columns of M is called proper if 

• None of the entries of M 2 is zero and 

• It contains exactly two distinct variables. 

Example 1: Consider the following matrix in three real 
variables xq,x\ and x 2 

Xo —Xl ~X2 

Xl Xo — X2 

X2 Xo Xl 

X2 —X\ Xo 



(5) 



The sub-matrix 



is proper while 







61 + 6 2 _= 61 © 6 2 © 6162, 
61 = 10 61. 



(4) 



X\ ~x 2 

X 2 X\ 

If M(i,j) ^ 0, then we" write \M(i,j)\ 

M(i,j) = Xk or — Xk ( in case of CODs, \M(i, j)\ 

M(i,j)€{±x k ,±x%}). 



is not. 



X3 

x 3 
= k whenever 



k if 



If M is a ROD and if M2 is a 2 x 2 proper sub-matrix 
of M , containing two variables, say xi and x m , I, m positive 
integers, then M% M% = [xj + x^)^- In other words, M-i is 
a ROD by itself in two variables. The following lemma gives 
a characterization of RODs in term of proper 2x2 matrices. 

Lemma 1: Let M be a p x n matrix in k real variables 
Xq,X\,X2,- ■■ ,%k-i, sucn that each non-zero entry of the 
matrix is Xi or — xi for some % € {0, 1, ■ • • , k — 1}. Then 
the following two statements are equivalent: 

1) M is a ROD. 

2) (i) Each variable appears exactly once along each column 
of M and atmost once along each row of M, 

(ii) if for some i,j,f, M(i,j) ^ and M(i,f) ^ 0, 
then there exists i' such that \M(i,j)\ — \M(i',j')\ and 
\M(i,f)\ = \M(if,j)\, 
(iii) any proper 2x2 sub-matrix of M is a ROD. 



A. STEP 1: Construction of a new class of square RODs 

Square RODs have been constructed by several authors, 
for example, Adams et al ||9) and Geramita et al ||T2l . All 
these designs are constructed recursively and the basic blocks 
of these designs are the RODs of order 1,2,4 and 8. It 
is known that these designs are obtained by left (or right) 
regular representations of the field of real numbers, the field 
of complex numbers, the Quaternion algebra and the Octonion 
algebra respectively. In this subsection, we take a different 
approach towards the construction of RODs and that lead to 
a new class of RODs constructable recursively of which the 
constructions of [9| and |12| are special cases. 

If B t is a square real design of size [t, t, k] in fc real variables 
xq, ■ ■ ■ , Xk-i, then whenever B t (i,j) 7^ 0, we write 



Theorem 3: Let t = 2 a and B t be a real design with 



(6) 



B t (i,j) = l^t(i,j)xx t (i,j), 
for some Ht{i,j) € {1,-1}, 

and \ t (i,j) £ {0,1,- ■■ ,k-l] 

for < i,j < t- 1. 

B t is uniquely determined by /i t and X t . 



The approach we take is identifying a pair of functions /it 
and At that defines a square ROD. Towards that end, for any 
t, identifying Z%a with FSji, we have S C Z^a identified with 
a subset of F§. We define two maps which are used for the 
construction of our square RODs as follows. 
Let 



It ■ Z p ( t ) -» Z t 



(7) 



be an injective map defined on Z p ^ with the image denoted 
by Z p(t) = lt (Z p(t) ) and 



If.Z, 



p(*) 



Zt 



(8) 



be another injective map defined on Z p ( t y 

In the following theorem, we define the two maps /i t and At 
given in (O which define the matrix B t , in terms of the two 
maps (O and ([H} and identify the conditions for the resulting 
B t to be a square ROD. 



B t (i,j) be non-zero if and only if i © j 6 Z 
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When 



't(i,j)^0, let 
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_ (_]\|*-Vt(*©j)l 


(9) 


>>t{i,j) = 


= 7t _1 (*©i)- 


(10) 



If 



\(ip t (x)®Mv))-(x®y)\ 



(ii) 



is odd, for all i,j/£ Z p / t \,x 7^ y, then, B t is a square ROD 
of size [t, t, p(t)}. 

Proof: We use the Lemma Q] to prove that B t is a ROD. 
First, for a fixed j £ Z t , define 

A = {i® j \i€Z t ,i®j€Z Kt) }. 

It is clear that A = Z„u\. Moreover, as 74 is injective, we 
have 7t~ 1 (i © j) 7^ l^ l {i' © j) whenever i ^ i'. Therefore, 
each column of the matrix B t contains all the variables 
Xq, Xi, ■ ■ ■ x p ( t )-i an d these variables appear exactly once. 
Similarly, it follows that the variables appear atmost once in 
any row of B t . 

Secondly, assume that B t (i,j) 7^ and B t (i,j') 7^ 0, then 
we show that there exists i' such that 

\B t (i,j)\ = \B t (i',j% 

\B t (i,j')\ = \B t (i',j)\. 

Let i' = i © j © j'. We have 

\B t {i, ] )\= 1 i 1 {i®i), 
\BS',3')\=lt 1 (i'® ] l )= 1 i 1 {i® ] ). 

Therefore, \B t (i,j)\ = \B t (i',j')\. Similarly, \B t (i,f)\ - 
\B t (i' ,j)\. Thirdly, we show that any proper 2x2 sub-matrix 
of B t is a ROD. It is enough to prove that /i t (i,j) ■ ^t(hj') • 
Vt(i',j) ■ IH{i',3') = -1' or equivalently, 

\i-Mi®j)\Mi-Mi®f)\ + \*'-Mi'®J)\ + \i'-Mi'®f)\ (12) 

is an odd number. But i' = i © j © j' and \x © y\ = |x| + \y\ 
where k = I if k — I is a multiple of 2. We can write (fT2l as 

\{i®i')-{tp t {i®3)®A{i'®3))\ 

= I ((*©i)©(*' ©j)) ■(&(»© i) ©&(»'© j)) I- 

which is an odd number as both i®j and i'®j are the elements 
of Z p (t). ■ 

Thus, it is enough to construct 7t and ijjt satisfying the 
property stated in Theorem [5J in order to construct a square 
ROD. This we achieve by making use of another set of two 
maps 4>i and fo as follows. 

On Z$ = {0, 1, ■ • • ,7} we define a map <f>\ : Z% 1— > Z& 
given by 



12 3 4 5 6 7 
12 3 4 7 5 6 

For a e {0, 1, 2, 3}, we have Z 2 « C Z 8 . Define 



02 



X I— > 01 (x) 



(13) 



(14) 



where the 2's complement of an element x of F2 is performed 
in Ff. Note that the map ijj 2 a is well defined even though it 
is defined in terms of </>i which is defined on Zg. 

Lemma 2: Let x, y £ Z 2 <^, a £ {0, 1, 2, 3}, x 7^ y. Identify 
Z 2 « with Ff. Then |(-02»(x) 8 ^2«(j/)) ■ (a: © |/)| is an odd 
integer. 

Proof: We prove it only for a = 3, For all other values 
of a, one can prove it similarly. Write x — ip$,{x). Let the 
radix-2 representation of x £ Z$ be (X2, X\, Xo) and that of x 
be (x2,ii,£o) where each Xi or Xj takes value from the set 
F2 = {0,1}. The following table describes the map x — ► x 
for x = 0, • • • , 7. 
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Using this table, and dU we express Xi in term of X2, X\ 
and xo for % = 0, 1 and 2 as follows. 

x — x 2 x + x 2 (xi © x ) — xo © X1X2, 

Xi = X 2 X\ + x 2 (xi © X ) — X © Xl X2X0, 

x 2 = S2a;i + xi(a>2 ©#0) = aro ©a>i ©X2 © xoaii . 
Hence 



,/. 1 



• 2_. x j © IT x ji i = 0, 1 and 2. 



j=o 



0<i<2 



Let x,y £ Zg, x ^ y. Now |(x © y) ■ (x © y)| is odd if and 

only if 

2 

y^(z» ©j/i)(£» ©&) = l. 
But Ei=o(^ © y<)(& © &) = l © nLo( 1 © ^ © yi)- 

Now 1 © IIi=o(l © Xi © y») i s equal to if an d on ly if 
1X7=0(1 © x » © y») = 1 i- e -' ( x * © y») = f° r a H *> which 
implies that x — y. 

As x ^ y, we have 5Z i=0 (x, © ?/*)(& © &) = 1. ■ 

The second map <fi 2 is defined on the set F given by 

F = ix £ V\ \x\ = l or 3} = {1, 2, 4, 7, 8, 11, 13, 14} (15) 

as injective map <f) 2 : F ^ Z w given by 



12 4 7 
12 4 6 



11 13 14 
15 10 12 



(16) 



Lemma 3: Let F be the set given in ( TTBT ) and 
x,y £ F,x 7^ y. Then 
(i) I < p2 (x) ■ x\ is odd for all x 7^ 0. 
(ii) |02(^) • y\ + \(j) 2 (y) ■ x\ is odd for all x ^ y, x ^ 0, y ^ 0. 



Proof: There are only finitely many possibilities for x and 
y and it can be checked that both the statements (i) and (ii) 
hold for all possible cases. ■ 

Now, we define the maps j t and ip t in terms of the maps 
0i and <f>2 as follows: The map 74 defined over Z p M is given 
by 



lt{i) 



where 



7(771) 



if < i < 7 
if i > 8, i = 8? 



m, < m < 7 



(17) 



7 = 



5 

11 



G 
13 



7 
14 



(18) 



For an element cc € ^p(t)> either x e Zg or 1 = 2 4y ~ 1 z for 
some j/6N\ {0} and z £ F. Let be the map defined on 
the set Z p {t) given by 



(x) = 



<j>i{x) 
2 4y " 1 • Mz) 



if x £ Z$ 
i£x = 2 4 y- 



l z, z£ F 



The map ipt is defined by 



M x ) = 0(«) in F^ Vx G Z p{t) 



(19) 



(20) 



The following theorem shows that the maps j t and ipt 
defined by ( fTTl ) and ( f20b satisfy the conditions of Theorem 
[3] and hence define a ROD. 

Theorem 4: Identify Z p (t) with a subset of Ff, £ = 2". Then 
|(^t(x)ffi Vt(y)) • (zffiy)j is odd for all x,y e Z p{t) ,x ^ y, 
and hence from Theorem [3] the matrix B t defined by j t and 
ipt by ([T7J and d20j is a square ROD. 

Proof: For £ = 1,2,4 and 8, the statement holds by 
Lemma |2] Hence we assume that t > 16. As ^t(0) = 0, 
it is enough to prove that 
(i) \ip t (y) ■ y\ is odd for all y ^ 0. 

(ii)IV't(ic) • y\ + \ipt(y) • x l is odd for a ii i t^i/ ! 

To prove (i), let z — ipt{y) ■ 2/- If y € ^> we have 
IV't(y) • 2/1 == (^(j/) • y| which is an odd number by Lemma 

rj 



On the other hand, if y 



)4Z-1„ 



I > 0,m £ F, then 



\z\ = |2 4/ ~ 1 02("i) • 2 4L ^ 1 m\ where the 2's complement of an 
element is performed in F|. We have \z\ = \(j)2(m) -m\ where 
the 2's complement of 4> 2 (m) is performed in F|. Hence \z\ 
is odd by Lemma [3] 

In order to prove the part (ii), we have following three cases: 
(i) 1 < X < 7 & 1 < y < 7, 

(ii) l<2/<7&2! = 2 ia - 1 f3 for some £ F, a > 1, 
(iii) x = 2 4 "- 1 /3 & y = 2 4q - 1 /3 for some ^, /3 £ F, a, a > 1. 
In all the three cases, we have x 7^ y. By Lemma|2] (i) is true. 
For the second case, let z — ip t (x) ■ y © ipt (y) ■ X. We have 
z = ( 2 4 "-y 2 (/3)' ■ y) © ((2 4 "- 1 /?) ■ 0T(y))- ' 

As 2 4a -V2(/3) -2/ = (the a ll zero vector in Ff) for a > 1, 
we have z = (2 4 "- 1 /?) • 0i(y). But |/3| is odd for all £ F, 
hence \z\ is an odd number. 

For (iii), let z — ipt(x) • y © ipt (y) ■ x. We have 



z = 2 4 «- 1 2 (/3) • 2 4& ~ 1 © 2 4 "- 1 /? • 2 4 «- 1 2 (/?)- 



If a > a, we have 2 4q ~ 1 /3-2 4 "- 1 02(/3) = and 2 4 «- 1 2 (/3)- 
2 4 "- 1 /3 = (3. Thus \z\ 
a = a, it follows that 



is an odd number by Lemma [3] If 



- 



\MP)-P\ + \P-MP)\ 



which is an odd number by Lemma [3] 

From Theorem [3] it follows that the matrix B t defined by 
7t and ipt by (fTTI i and ( f20b is a square ROD. ■ 

The square RODs of Theorem [4] will be denoted by Rt 
throughout. The ROD R\q of size [16, 16, 9] is given by (|2"H 
at the top of the next page. As another example the ROD R32 
of size [32, 32, 10] is given by d22l in the following page. In 
Appendix [A] it is shown that the RODs R t can be constructed 
recursively. 

One can define the functions -f t and i/j t different from the 
one given above and can have a square ROD different from 
R t . In Appendix iBl we provide three different pairs of such 
functions and these are shown to give the well known Adams- 
Lax-Phillips' construction from Octonions and Quaternions 
and Geramita and Pullman's construction of square RODs. 

B. STEP 2 : Construction of new sets of rate-1 RODs 

Transition from a square ROD to rate-1 ROD is illustrated in 
J7) using column vector representation of a ROD. In a similar 
way, we construct a rate-1 ROD W n of size [v(n),n, v(nj\ for 
n transmit antennas from a ROD of size \v{n), v(n), n] where 
n is any non-zero positive integer, not necessarily power of 2. 

Any square ROD of order v(n) obtained via a suitable 
pair of mapping j v ( n ) and ip v t n ) satisfying the conditions 
of Theorem [3] (for instance, i?„( n ) obtained in the previous 
subsection or A v ^ n ), A u i n \ and G v r n \ obtained in Appendix 
IB1 can be used for this purpose. We refer to any such design 
by 5„( n ) consisting of n real variables Zq,Zi, ■ ■ ■ , 2,1-1. 

Let yo,yi, ■ ■ ■ ,y u ( n )-i be v{n) real variables which con- 
stitute the matrix W n . The matrix W n is obtained as follows: 
Make W n (i, j) = if the i-th row of -B„( n ) does not contain 
Zj. Otherwise, W n (i,j) = y k or -y k if B„(„)(i, k) = Zj or 
— Zj respectively. The construction of the matrix W n ensures 
that it is a rate-1 ROD. Using ©,(|9]) and Theorem[4] we have 



where 



W n (i,j) = 8(i,j)yf(i d ), 
s(i,j) = (_i)W«c»)('M»)tf)}l 

f(i,j) = i®7v(n)(j) 

for < i < v(n) - 1, < j < n - 1. 

Let W n be a matrix which is similar to the matrix W n 

by 



W„(i,j) = 8(i,j)yf(i d ), 



(23) 



(24) 



given 



(25) 



where 



s(i,j) = (-ljl^wW^-w^wW)! (26) 

f{t,j) = i®7v{n)(j)- 

W n is also a rate-1 ROD. Both W n and W n are the new sets 
of rate-1 RODs that are used in the following subsection to 
construct our codes. 



As examples, the RODs Wg and Wg of size [16, 9, 16] 
obtained using it!i6 are given by d27l i and (f28b respectively 
and the RODs W w and W w of size [32, 10, 32] obtained using 
R32 are given by J29] i and (T30b respectively. 
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J/1 -J/14 -J/22 
J/0 -J/15 -J/23 
J/15 J/0 -J/24 
J/1 -J/25 
J/2 -J/26 
J/3 -J/27 
J/4 -J/28 
J/5 -J/29 
J/6 -J/30 
J/7 -J/31 
J/23 J/24 J/0 
J/22 J/25 J/1 
J/21 J/26 
J/20 J/27 
J/19 J/28 
J/18 J/29 
J/17 J/30 
J/16 J/31 
J/31 -J/16 
J/30 -J/1 7 
J/29 -J/18 J/10 
2/28-2/19 J/11 
2/27-2/20 J/12 
J/26 -J/21 J/13 
J/25 -J/22 J/14 
J/24-J/23 J/15 



J/2 
J/3 

J/4 
J/5 
J/6 
J/7 

J/8 
J/9 



(30) 



C. STEP 3 : Construction of low-delay rate- 1/2 Scaled CODs 

In this subsection, we construct a rate- 1/2 scaled COD 
with the help of rate-1 RODs W n and W„ constructed in the 
previous subsection. 

Let Xq,Xi,--- be complex variables. The 8x8, rate-i 
CODs A(xo,Xi,X2,X3) and B(xi,xs,xe,x 7 ) and the 8 x 1 
column vector C(xo,xi,X2,X3) in four variables, shown be- 
low, are the basic ingredients for our construction of rate-i 
scaled CODs. 



A(xo,Xl,Xl,X3) 



:/■()- 

Xl 
X2 



X3 







x, 



xi- 


X3 






■' n 



B(x4,X5,x 6 ,x 7 ) 



COo,Zl,IE2,Z3) = -7= 
V2 



X4- 


-x,- 


-x fi - 
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'''4 
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-Xs 
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''4 





X 7 
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-x 5 








X 7 - 


-X6 















•"1 
xo 

J(l 

0-IEl 

0-x 2 

X3 0- 





xi 








4 

0-^3 

0- 

c l x 2 


x - 

XI 



xo 


X2 



0- 



x 5 

X4 





0- 




'■3 



'% 

r\ 
'■0- 




x 7 

0— Xj 
xl 

x* 
X4— xl 
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X7-X6 x: 5 x 4 



X01I0 lU-i J>Qj>r\ tO^ 3^2 **s 3 



(32) 



(33) 



Let i be any non-negative integer. Define the following 
matrices: 

A(2i) = A(x$i,x$i + i,x 8 i+2,xsi +3 ), 
A(2i + 1) = B(x Sl+ 4,x Sl+5 ,x 8 i +6 ,xs l+7 ), (34) 

One can easily verify that the matrices given by 



A(0) 
A(l) 



3(1) 
A(Q) 



A(l) 
A(Q) 



-A(0) 
A(l) 



(35) 



are scaled CODs as the columns of the matrices are orthogonal 
to each other and the norm of each column is equal to square 
root of the sum of the norms of the variables of the design. 
By relabeling the variables in the matrices A(0) , A(l) , A(0) 
and A(l), it follows that 



A(i) 



A(i) 



is a scaled COD whenever (i + j) is odd and 



A(i) 

A{j) 



A(i) 



(36) 



(37) 



is a scaled COD for all values of i and j, i ^ j. 

Let n be an integer such that n > 9. We construct a matrix 
(DR) n of size v(n) x n as follows: Let t = n — 8 and 
Wt and Wt are the two rate-1 RODs of size [i/(t),t, v(t)] 
in v(t) real variables yo,yi,--- ,y v (t)-i constructed in the 
previous subsection. Let H t and H t be the matrices formed 
by substituting y, with A{2i + 1) in the matrix W t and A(2i) 
in the matrix Wt respectively for i = to v{t) — 1. Note that 
the size of both H t and H t is 8v(t) x t. 

Let u = v{n)/8. Let E% and Og, each of size 4u x 8, be 
defined as follows: 



#8 = 



A(0) 
A{2) 



o 8 = 



A{u-2) 
Define the matrix (DR) n as 

(DR) n = 



Eh 
Oh 



A(l) 
A(3) 



A(u-l) 



H t 



(38) 



(39) 



Note that the number of rows and columns of the matrix 
(DR) n are 16 • v(n - 8) = 8 • v(n)/8 = ^(n) and i + 8 = n 
respectively. The following theorem is the main result of this 
paper. 

Theorem 5: Let n be a positive integer and (DR) n be the 
matrix as defined in d39l . Then (DR) n is a rate-1/2 scaled 
COD of size [i/(n),n,^]. 

Proof: For n < 8, one can construct rate-1/2 COD of size 
[v(n),n, 4^] from a COD of size [8, 8, 4] given in (ED- We 
assume that n > 9. Let p = f(n). We have 



(DR)%(DR) n = 






E s 



HfEs + HfO; 



ofo 8 



E s Ht + O s Ht 
H t H H t + H t H H t 



From the construction of Eg, and O s given in ( 1381 ), we have 
E*Ek + OgO s = (|x | 2 + • • • + Ixv^h. 
From equation d37| i, we have 



H?H t + H t H H t = (|xo| + ' ' ' + k P /2-i| )4- 



Thus it is enough to show that E% H t + 0% H t = Og X ( n _g) 
where 0$ x (n-8) is a matrix of size 8 x (n — 8) containing 
zero only. Let the j'-th column of H t and H t be H t (j) and 
Ht(j) respectively. Then we show that Z(j) = E% H t (j) + 
0£H t (j) = 8xl for all j e {0, 1, • • • ,n - 8 - 1}. 
Let u = p/8. For convenience, we write 7 for 7„(t). We have 



E» = 
0§ = 



A H (0) A H {2) 
A H {1) A H (3) 



A H (u 
A H {u 



Ht(j) 



S (0,j)A(2(0e 7 (j)) + l) 
s{l,j)A(2(l © 7 (i)) + l) 



s(z,j) J 4(2(i©7(i)) + l) 



9 (|_l,j)^(2((|-l)© 7 (i))+l) 

s(0,i)3(2(0© 7 (i))) 

s(l,j)4(2(le7(j'))) 



*(m)A(2(i e 7O'))) 



|(|-l,j) j4 (2((|-l)©7(i))) 



where s(i,j) and s(i,j) are defined in d23b and (1251 ) respec- 
tively. We have 



fftCi) 



^(i) = ^^,i)A ff (2^(2(i© 7 (i)) + l) 

i=0 

f-i 
+ £ s(i,i)A H (2i + l)A(2(i © 7 (i))). 



Now 



X)S(*,i)^(2t+l)A(2(t© 7 (3))) 

= ^ s(i © 7C?),i)^ H (2(t © 70')) + 1)^(2*)) 
and s(i,j) = s(i®j(j),j). 



i=o 



Hence 

Z(j) 



£(*(*,i)^(2t)A(2(t© 7 (j)) + l) 

i=0 

+§(* © l{]),])A H {2(i © 70')) + W) 
f-i _ 

53*(*,i)(^(2i)A(2(t©7(j)) + l) 



j=0 



+A h (2(i®j(j)) + 1)A(2i)) 
8x i using 



We illustrate our main result in the following example. 
Example 2: For 9 transmit antennas, we have rate-1/2 
scaled COD of size [16, 9, 8] given by 
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v'l 

~- T '.-l 

■'•"1 
A 

~ x 3 

vf 

;r 2 
v" 
~ x l 
j/2 

^f" 

' V) 

7 

- x l 

v5 

'2 



x O 
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while the known rate 1/2 scaled COD for 9 transmit antenna 
is given by @ 



->- 
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where the decoding delay is 32. 

For 10 transmit antennas, the rate-1/2 scaled COD given by 
Tarokh et al [7| of size [64, 10,32] is given in Appendix ICl 
while the new rate-1/2 code of size [32, 10, 16] is given by 

(ED. 
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(41) 



D. Maximal rate of the scaled CODs 

It has been shown by Liang [ 1 1 that the maximal rate of a 
COD for n transmit antennas is i + ^ when n = 2t — 1 or 
2t 

The following result says when there is scaling of atleast 
one column then the maximal rate is i. 

Theorem 6: The maximal rate of a scaled COD, with scal- 
ing of at least one column, for n transmit antennas is \. 

Proof: Let (DR) n be a scaled COD for n transmit 
antennas and there exists at-least one column of the matrix 
such that whenever a variable appears in that column, it is 
scaled by -k= . Since all the variables appearing in a column is 
scaled by -h=, each variable must appear twice in that column. 
Let the number of distinct complex variables in (DR) n is k. 
Then 2fc < p, i.e., k/p < 1/2. ■ 

Thus if one allows to incorporate a factor of -4= for all 
the entries of any column, there won't be any improvement in 
the rate, on the other hand, as we will have observed, we can 
construct some rate 1/2 codes with lesser decoding delay. 
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III. Delay-minimality for 9 transmit antennas 

In this section, it is shown that the low-delay rate 1/2 COD 
developed in the previous section is of minimal delay for 9 
transmit antennas. To prove this, we need some preliminary 
facts regarding the interrelationship between real and complex 
ODs and certain bilinear maps. The ROD and bilinear maps 
are intimately related in the sense that a ROD of size [p, n, k] 
exists if and only if there exists a normed bilinear map of 
the same size. The normed bilinear maps have been studied 
extensively and one find a good introduction to this topic in 
the book by Shapiro ifTSl . As the results from the theory of 
normed bilinear maps is used to prove our claim, these maps 
have been defined below and some facts are stated regarding 
these maps. 

A bilinear map / (over a field F) is a map 

/ : ¥ k x F" -» ¥ p (42) 

(x,y) ^ f(x,y) (43) 

such that it is linear in both x and y, i.e., f(xi + X2,y) = 
f(xi,y)+f(x 2 , y) and f(x, yi +y 2 ) = f(x, yi)+f(x, y 2 ) for 
all x,xi,x 2 <E F fe and y,yi,y2 G F n . The space F p is called 
the target space of /. If the vector spaces under consideration 
are inner product spaces, for example, when the field is real 
numbers or complex numbers, the Euclidean norm of a vector 
X is denoted by || x ||. If a bilinear map preserves the norm, 
then it is called a normed bilinear map. More precisely, 

Definition 3: A normed real bilinear map (NRBM) of size 
[p, n, k] is a map / : R fe x R™ — j- W such that / is bilinear 
and normed i.e., || f(x,y) \\ = \\ x \\\\ y | Vx <E R k ,y € R". 
If f(x,y) = implies x = or y = 0, then such a map is 
called a nonsingular map. 

The following theorem gives a lower bound on p for fixed 
values of n and k. 

Theorem 7 (Hopf-Stiefel Theorem [18]): If there exists a 
nonsingular bilinear map of size [p, n, k] over R, then (x + 
y) v = in the ring W 2 [x,y]/(x n ,y k ). 

Definition 4: Let n, k be positive integers. Then the three 
quantities no k, pbl and pnbl are defined by 

• no k — min{p : (x + y) p — in F2[x, y]/(x n ,y k )}, 

• PBh(n, k) — min{p : there is a nonsingular bilinear map 
[p, n, k] over R }, 

• PNBL{n, k) ~ min{p : there is a normed bilinear map 
[p, n, k] over R}, 

The following basic facts about these quantities are well 
known fT8l . 

Pnbl (n, k) > p B L {n, k) > n o k and pnbl (n, k) = n if and 
only if k < p(n) where p is the Hurwitz-Radon function. 

It follows from the definition of n o k that 

Proposition 1 ([18]): no k is a commutative binary opera- 
tion. 
(/) If k < I then no k < no I 

(II) n o k = 2 m if and only if k, n < 2 m and k + n > 2 m . 

(III) If n < 2 m then no(k + 2 m ) = nok + 2 m . 
Example 3: Let us compute 10 o 10. We observe that 10 < 

2 4 , but (10 + 10) > 16. So, 10 o 10 = 16. 

The relation between RODs and NRBMs has been observed 
by Wang and Xia in |17|. For the sake of completeness and 



since the proof of this fact gives the explicit relation between 
the NRBM and the row-vector representation matrices of the 
ROD which is in correspondence with it, we give here the 
proof of the following lemma. 

Lemma 4: A ROD of size \p, n, k] exists if and only if there 
exists a normed real bilinear map of size [p, n, k]. 

Proof: Let x = (x\, x 2 , ..., x k ) T , y = (yi,- ■■ ,y n ) T and 
z = (-?!,•■• , , z p ) T be real column vectors. 

Let A be a ROD of size [p, n, k] given by A = J2i=i ^i x i 
where the p x n real matrices are the dispersion matrices or 
weight matrices defining the design A QJ. Let 



/ :R fe x 



{x,y) i-> (^2AiXi)y. 



We show that / is a normed real bilinear map of size [p, n, k\. 

Let z = f(x,y). Then z% — xP Biy where the k x n real 
matrices Bi, i = 1, 2, • • • ,p, are the row vector representation 
0] of the design A. This representation shows that / is 
bilinear. Moreover, / is normed, since || f(x, y) \\ 2 —\\ Ay || 2 = 
(Ay) r Ay= yj (x\ + x\ + ... + 4)I n )y =|| x l| 2 || V \\ 2 ~ 

To show that the converse holds, let / be the normed bilinear 
map given by 



/ 



X M" 

(x,y) 



As / is linear in both x and y, we have z = Ay where A is 
anpxn matrix where each entry of the matrix is a real linear 
combination of the variables xi, • • • ,Xk- As / is normed, we 
have || z || 2 =|| f(x,y) || 2 =|| x || 2 || y || 2 . But f(x,y) = Ay. 
Then, || Ay || 2 = ((x\ + ■ • • + x\)I n )y T y. So, we have 

y T Ay =((x\ + ■■■ + -. xt)I n )y T y. 

As y consists of variables, the above equation is equivalent to 
A T A= (x\ +xl + ... + x 2 k )I n . ■ 

We now prove the main result of this section. 

Theorem 8: The minimum decoding delay of the rate-1/2 
COD admitting linear combination of complex variables for 9 
transmit antennas is 16. 

Proof: In this proof, we assume that the ROD and COD 
admit linear combination of two or more variables as its 
entries. Suppose, there exists a COD of size [2x, 9, x] where 
x is an integer less than 8. This implies existence of a ROD 
of size [4x, 18, 2x},x < 8 which is obtained by replacing each 
complex entry by its 2 x 2 real matrix representation. In the 
remaining part of the proof, we show that such a ROD does 
not exist thus proving the theorem. 

If a ROD of size [Ax, 18, 2x] exists, then there also exists a 
normed bilinear map of the same size by Lemma [4] which 
implies that 4a; > 18 o 2x. As 18 o 2x > 18, we have 4x > 18 
i.e., x > 5 for x being an integer. Thus we have three possible 
choices for x, namely 5, 6 and 7. 

By Proposition[U we have 18o2a; = 26, 28, 30 for x — 5, 6 and 
7 respectively. In all the cases, 18o2ir > Ax which contradicts 
the fact that Ax > 18 o 2x. ■ 
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IV. PAPR REDUCTION OF RATE- 1/2 SCALED CODS 

In this section, we study the PAPR properties of Scaled 
CODs constructed in this paper. Note that in the construction 
of TJC n given in Q, even though the delay is more, there are 
no zero entries in the design matrix. On the contrary, in our 
construction of low-day codes (DR) n there are zero entries. 
To be specific, observe that the first eight columns of rate-1/2 
code (DR) n ,n > 9 given in (|39l contains as many zero as 
the number of non-zero entries in it, while there is no zero 
in the remaining columns of the matrix. When the number 
of transmit antennas n is more than 7, the total number of 
zeros in the codeword matrix is equal to 8(v(n)/2) = 4i/(n). 
Hence the fraction of zeros in the codeword matrix is equal 
t M=) = 4/n for n > 8. 

Now in the remaining part of this section, we show that one 
can further reduce the number of zeros in (DR) n by suitably 
choosing a post-multiplication matrix to it without increasing 
signaling complexity [15| of the code. 

As seen easily, only the first eight column contain zeros 
while the others do not. Moreover, the zeros in 0-th column 
and the 7— th column occupy complementary locations, so is 
also for the pairs of columns given by (1, 6), (2, 5) and (3, 4). 
What it essentially suggests is that we can perform some 
elementary column operations which will result in a code in 
which all the entries are non-zero. In other words, if the rate- 
1/2 COD is of size p x n, then we post-multiply it with a 
matrix Q n of size n x n given by 

A 
I n -i 



^cn — 



where A is a matrix of size 8x8 given by 



V2 



and I„_8 is the (n — 8) x (n — 8) identity matrix, then all the 
entries of the scaled COD given by (DR) n Q n , are non-zero. 
We formally present this fact as: 

Theorem 9: The matrix Q n when post-multiplied with a 
rate-1/2 scaled COD (DR) n given by 
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Eg 
O* 



H t 



(44) 



always gives a COD with no zeros. Moreover, the matrix 
Q n does not depend on any particular construction procedure 
(namely the maps 7 t and ip t ) used to obtain the constituent 
rate-1 RODs. 

Proof: It is clear that the first 8 columns of the matrix has 
50% zeros in it and in the remaining n — 8 columns formed 
by H t and H t , there are no zeros as both these matrices are 
constructed from rate-1 ROD by substituting all the variables 
in it with appropriate 8-tuple column vectors. Here neither 
rate-1 ROD nor the 8-tuple column vector has any any zero 



in it. Therefore, the matrix Q n gives a rate 1/2 scales COD 
without any zeros irrespective of how the rate-1 RODs are 
obtained for the construction of (DR) n . ■ 

Example 4: The rate-1/2 code with no zero entry for 9 
transmit antennas is given by 
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-2/5 


2/4 


-2/4 


-2/5 


2/6 
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2/5 


ill 


.'4 


-2/1 


-2/6 


2/7 
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2/7 
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-2/5 
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-2/5 


-2/1 
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2/6 


-J/5 


-2/4 


-J/3 



V. Discussion 

This paper gives rate-1/2 CODs for n transmit antennas 
with decoding delay equal to v(n). The decoding delay of 
these codes is half of that of the rate-1/2 CODs given in Tarokh 
et al 0. As the maximal rate of a scaled COD is close to 1/2 
for large number of transmit antennas, the codes constructed 
in this paper are better than the codes constructed by Liang [ 1 ] 
or Lu et al ifTUI when considered for large number of transmit 
antennas. Another advantage with the designs reported in this 
paper is that they do not contain zero entries leading to low 
PAPR. 

All the four constructions namely Adams, Lax and Phillips 
construction from Quaternions, Octonion, Geramita-Pullman 
construction and the construction given in this paper will give 
the same square ROD if number of transmit antennas is less 
than or equal to 8. Therefore, these four construction will 
generate the same rate 1/2 scaled COD if the number of 
transmit antennas ( of the scaled COD) is less than or equal 
to 16. For more than 16 antennas, rate-1/2 scaled CODs will 
vary with the methods chosen for the construction of rate-1 
RODs. Due to the largeness of the matrices involved, it is 
not possible to display two distinct rate-1/2 scaled CODs for 
17 transmit antennas, obtained by two different construction 
procedures for rate-1 RODs. 

It is not known whether the low-delay for rate 1/2 scaled 
CODs we have achieved is minimal delay except for the 
case of 9 transmit antennas. We conjecture that v(n) is the 
minimum value of the decoding delay of rate-1/2 scaled CODs 
for any n transmit antennas. It will be interesting to see 
whether this is indeed true. 

An interesting direction for further research would be to 
investigate whether the necessary conditions given in Theorem 
[3] on the maps j t and ipt are indeed sufficient also. 
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It follows that 
Kj = K?(x ,x 1 ,- 



and 



Kl 



,xt-i) = K t (x ar 
= K t (-x ,Xi,-" 



,X t -l) 



for t = 1, 2,4 or 8. The expression for Rt of order t as given 
in Theorem|4]gives rise to the following recursive construction 
of R t . Given two matrices U — (utj) of size v\ x w\ and V 
of size V2 x W2, we define the Kronecker product or tensor 
product of U and V as the following V1V2 x W1W2 matrix: 

I U11V ui 2 V ■■■ u 1Wl V ^ 
uuV U12V ■■■ u lw ,V 



\ u Vll V u Vl2 V ■■■ u VlWl V ) 
Let I n be an identity matrix of size n. Define 
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h, 
II 



In ® I 
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be real variables. Define 



Vili 



ysil- 



Let 2/0, • • • ,y 5 

24(2/0, yi) = yoi2 + yiiL 
T 8 (y2, 2/3, 2/4, 2/5) = y^lg + yalg - 

We have four RODs of order n = 2 a with a = 0, 1, 2, 3 as 
given in (l45l l which are respectively K\, K2, K4 and .Kg. 
Assuming that a square ROD of order n = 2 4l ~ 1 , 1 > 1 

Rn = Rn(xo, ' ' ' , X p (n)-l) 

which has p(n) real variables, is given, then we construct 
R2m Rim Rs,m Ri6n of order In, An, 8n and 16n respectively 
given by d46l i. as shown at the top of the next page where 



Vi — x p(n)+2+i 

Rt 



and 



T 



-rT 



RT(xq,xi,--- 
R t (x ,-xi,-- ■ 
Rt(-xo,Xi,-- ■ 



x p(t)-i) 
:— x p(t)-i), 

: X p(t)-l)- 



Appendix A 
Recursive Construction of R t 

In this appendix we show that the RODs R t can be 
constructed recursively. 

Let K t = B t for t = 1, 2,4 and 8. The four square ODs 
Kt, t = 1, 2,4, 8 are shown below. 
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Appendix B 

Adams-Lax-Phillips and Geramita-Pullman 

constructions as special cases 

In this appendix we show that the well known constructions 
of square RODs by Adams-Lax-Phillips using Octonions and 
Quaternions as well as the construction by Geramita and 
Pullman are nothing but our construction corresponding to 
specific choices of the functions 74 and ipt defined by ^} and 
©. It turns out to be convenient to use the map \t = tytlt 
instead of the map ip t . Note that both j t and \t act on tne 
set Z p (j^ and are injective. Now given 74 and \t, we have 
ipt = Xtlt ■ With this new definition, we can reformulate 
the criteria given in Theorem |4] as follows. 



\{xt(x)®xt{y))-(it(x)®jt(y))\ 

is an odd integer Vx, y £ Z p (t)> x ¥" V- 



(47) 



13 



R2n = 



R$,n — 



R 



4)1 



T4,(—yo,yi)®I n 



**"n X p(n)J-n 

~Xp[ n )^-n J^n 

Ti{ya,yi) <^ I n 

I ^An 



Rin — 



Rl&n — 



Rln 

~ x p{n) + lI L 2n 



T s (-y2,ya,y4,ys) 



p(n) + l-'2n 



R 



T&(y2,yz,y4,ys) ® I n 



In R$n 



(46) 



In the following lemma, we define j t and \t m three different 
ways and these maps are shown to satisfy the relation given 
in d47b . Although both j t and Xt are different for all the three 
cases for arbitrary values of t, 74 is the identity map when 
t = 1, 2, 4 or 8. Hence \t = ipt if t G {1, 2, 4, 8} and is given 
by CH), 

Lemma 5: Let £ = 2 a , a = 4c + d, m G {0, 1, • • • , 7}. Let 
7t and xt be two maps defined over Z p ^ in three different 
ways as given below. Identify j t (Z p ^) and Xt{Z p (t)) as 
subsets of F§. Then |(7t(zi)®7t(z 2 )) • (xtOOeXtO^))! is 
odd for all Xi, X2 G ^(t), £1 7^ £2- 
For a; = 8/ + m G ^ p (t) , 

(i) 



7t(8/ + m) 

Xt(8l + m) 

(ii) 

7t(8/ + m) = 
X t(Sl + m) = 



where X4 = 
(hi) 

7t(8Z + m) 



Xt(8l + m) 
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; ) + 2 2 
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if Z 7^ c, m 7^ 0, 
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L ) + 2 i '(m-4) if 4 < m< 7, 



-2; 
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'1) 
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ifl^c,me {1,2,3} 
if I ^ c, m e {5, 6, 7} 



12 3 
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-4/ 
-4Z 







2 Z 




X2«( m ) 





m 



2 L 



tXsjrn) 

2 4(! + l) 



if I < c, 
if l = c 

if Z = 0,to = 
if Z ^ 0, to = 

if Z = c, m 7^ 0. 

if Z 7^ C, 777, 7^ 0. 



Proof: We give proof only for the case (i). The cases (ii) 
and (iii) can be proved similarly. 

It is enough to prove that 
(Bl) |7t(x) • Xt{x)\ is odd for all x 7^ 0, x G Z p (t) an d 
(B2) I'ftixi) ■ Xt(x2)\ + \lt(x 2 ) ■ Xt(xi)\ is odd for all xi,^ G 
Z p (t), Xi 7^ X2,X\ 7^ 0, x 2 7^ 0. 

Let 7t(8Z + 777) = 7 t (1) (8Z + m) + 7 t (2) (8Z + m) such that 



Similarly, let xt(8Z + 777) = x^ {81 + m) + x^ {81 + m) such 
that 



Xt (1) (8Z + m) 



P) 



X ( t'{8l- 




if Z = 0,777 = 0, 
if Z 7^ 0, 777 = 0, 

if Z = c, 777 7^ 0, 
if Z 7^ c, 777 7^ 0, 

if Z = 0, 777 = 0, 
if Z 7^ 0, m = 0, 
if Z = c, 777 7^ 0, 
if Z 7^ c, m 7^ 0. 



Let 8Z + 777 7^ and 8Z' + m! 7^ 0. From the definition of 
7J, xl' i = 1; 2, it follows that 

(Al) \x? ) {8l + 777) • 7t (2) (8Z' + m') I = if Z ^ V , 
{A2) | X «(8Z + 777) - 7t (1) (8Z' + m')l = 1 if I < I', 

{A3) |x| 1) (8Z + 777)-7 t (1) (8Z' + 777 , )l = if I > V 
or if I — I', m 7^ 0, 

(A4) Ix^ (8Z) - 7t (1) (8Z + m)| = l if Z 76 0, 

(A5) Ix^W-Tf^HlxfW-^^HO 

V x,y G Z p(t) , 

(A6) | X f } (8Z) • 7t (2) (8/ + ro)| = |x| 2) (8Z + m) • 7t (2) (8Z)| = 0. 
First we prove (Bl). Let x = 81 + m with 777 7^ 0. We have 

|Xt(aO -7*0*01 = |x^ 1) (8' + ™)-7 t (1) (8/ + m)| + |xf ) (8/ + m) 
■ 7t (2) (8/ + m)| + Ixt 1 ^ 8 ' + m ) ' l ( t 2) {81 + m)\ 
+|Xt 2) (8/ + m)- 7t (1) (8/ + m)| 
= | x ( 1 )(8« + m)- 7t (1) (8/ + m)| 

+|Xt 2) (8Z + m)-7 t (2) (8/ + m)| by (A5) 

= \\t ) (81 + rn) ■ 7 t (2) (8/ + m)\ using (A3) 
= |Xe(«7) ■ m\, e = 2 d if I = c, else e = 8 

But |x e (TO) • 777 1 is an odd number by Lemma [2] 
If m = 0, we have \-ft{x) ■ Xt{x)\ — 1 by (A4). 

To prove (B2), let x\ 7^ and x 2 7^ 0. Write X2 = 8Z2 +777.2, 
X\ = 8Z1 + 777i with X2 > x\. We have two cases: 
(CI): l 2 > h, (C2): l 2 = h = I, m 2 > mi. 

Case (CI): we have 



Xt(a>2) ■ 7t(»i) = Xt 1) (8«2 + m s ) ■ 7 t (1) (8h +mi) 
ext 2) (8« 2 +m 2 ) -7 t (2) (8«i +mi) by (A5) . 



.(1)/ 



J 1 ) 



7 t (1) (8Z 



i(l -2-') and7 t (2) (8Z 



8*777. 



But |xJ iJ (8Z 2 + 777 2 )-7r+8Zi+777i) 
and |x t (2) (8Z 2 + 7772) ■ if ] {8h + m x )\ 
thus \xt{x2) -lt{xi)\ = 0. 



= by (A3) 
by (Al), 
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Now xt(xi) ■ jt(x 2 ) = X t (1) (8^i + mi) • 7t (1) (8?2 + m 2 ) © 

xi a) (8li + mi) • 7t (2) (8^ 2 + ma) by (A5). 

But |x| 2) (8?i + mi) • 7 t (2) (8^ 2 + m 2 )| = by (Al) and 

\x?\Ui + mi) • 7t (1) (8^2 + m 2 )| = 1 by (A2). 

Hence |xt(»i) ■ Jt(x 2 )\ + \xt{x 2 ) • 7t(a?i)| is an odd number. 

Case (C2): we consider two following cases: 

(i) mi 7^ and (ii) mi = 0. Note that 7712 is always non-zero. 

Let d = \(xt(xi) • 7t(ar 2 )) © (xt(aj 2 ) • T*(«i))l- 

Case (i): We have 

d = |Xt 2) (8i + mi)-7 t (2) (8« + m 2 )| 

+|xf \%l + m a ) ■ li 2) (8l + mx)\ by (A3) and (A5) 
= l(Xe("ll) ■ ma) © (Xe( m 2) ■ mi)\, e = 2 d if / = c, else e = 8 

which is an odd number by Lemma [2] 
Case (ii): Since mi = 0, therefore I ^ 0. We have 

d = \xi 1) (8l)-li 2) (8l + m 2 )\ 

+ | x r ) (8/ + m 2 ). 7t (1) (80|by(A6). 
= 1 by (A3) and (A4). 

■ 
By Lemma [5] and Theorem [3] the matrix B t defined by the 
two functions 7t and xt is a square ROD in all the three 
cases. We refer to these three different RODs by A t ,A t and 
Pt corresponding to the pair of functions defined in (i), (ii) and 
(iii) respectively. Observe that our construction R t is different 
from any of A t , A t and P t for general values of t. 

Now, we proceed to show that the designs A t , A t and P t are 
essentially the Adams-Lax-Phillips construction using Octo- 
nions and Quaternions and the Geramita-Pullman construction 
respectively with change in sign of some rows or columns. 



Lemma [6] that A t is same as 



»(0) 



with t — 16n. 



Lemma 6: Let t > 16 and a power of 2. Also, let A t be 
the square ROD of order t as given in Lemma[5](i), and 16rl 
be the square ROD given in d49b which is of order \6n. Then 



A f = 



16n 



for £ = IQn. 



Proof: We prove it by induction on t. For t = 1,2,4 and 



A f = K t and the COD 



»(0) 



of order i is also given by 



if t . Hence the lemma holds for t = 1,2,4 and 8. Assuming 



that the lemma holds for t = n, i.e., A n = 



_<n»(°) 



of order n, 



we have to prove that the lemma also holds for t = 16n, i.e., 

16n 



AlSn 

Let 



.4 



16n 



An 
i 2 i 



A12 

A 22 



(50) 



where Ay, 1 < i, j < 2 are square matrices of size 8n x 8n. 
It is easy to check that the location of non-zero variables in 
the matrix Aiq h coincide with that of 0^ 6 „. Therefore it is 
enough to show the signs (positive/negative polarity) of the 
corresponding entries in the two designs are same i.e., 

1) Mi6n(*,i) = Mi6n(i%8,j%8) for < i,j < 8n- 1, 

2) /Xi6n(*,i) = Ms(i)i) for < i, j < 7, 



3) Mi6n(i,j) = Mi6 n (»®«%8,j© j%8) 
if < i < 8n - 1, 8n < j < 16n - 1, 

4) |Ui6„(8f, 8n © 8j) = Mn(«,i) for < i,j < n- 1, 

5) Mi6n(8n © i,8n © j) = Mi6n(«,i) if « © j = or 
« © j > 8n, 



A. Adams-Lax-Phillips Construction from Octonions as a spe- 
cial case 

The Adams-Lax-Phillips construction from Octonions is 
given by induction from order n = 2 a to 16n as follows QJ: 
Denoting the square ROD of order n — 2 a resulting from the 
Adams-Lax-Phillips construction using Octonions by 

O n = O n (*0,- •• ,Xp(n)-l) 

which has p(n) real variables, the square ROD of order 16n 
with (p(n) + 8) real variables x i: i = 0, 1, • • • , p{n) + 7, 



nen 



"w 



n(XQ, ■ 



o(n)+7) 



is given by 



I n (g) Kg(y , ■ ■ ■ , y-i) <Q n ®h 

®Z®h I n ®(-K[(y 0l - 



,2/7)) 



with y t = x p ( n)+i . 

With re-arrangement of variables and change in signs, we 

rewrite the design ©i6n as 



6) /ii6n(8nffii,8nffi j) = -/ii 6n (i,j) 
if i@3 e{l,2,-.- ,7}U{8n}. 



Note that 
1) & 2) together imply An = I n <g> K 8 (xo, ■ ■ • , xr), 
3) & 4) together imply A 12 = ol 0) <g> I s and 
5) & 6) together imply A 22 = A[ 1 ,A 2 i = -Aj 2 . 
Let Ai 6n (i,j) t^O. 

Then i® j e Z p(16n) and fnUhj) = (-1)^^®^ . 
To prove 1), we have to show that \i ■ tfji6 n (i (Bj)\ = |(i%8) • 
ipi6n(i%8 © j%8) for < i,j < 8n - 1. 
We have i®j= (16n)(l - 2~ l ) + 8 l m and i © j < 8n. So 
I = and i © j = m. i.e., i®j = i%8 © j%8. 
Thus it is enough to prove that \(i © i%8) ■ tpie n (i ® j)\ = 
Now (i © i%8) < 8n, 8 divides (i © i%8) and ^i6«(i © j) = 
8n © ips(m), hence the statement holds. 

The statement 2) is true as \i ■ ipi6 n (i ®j)\ — \i' ^s(* ®j)\ 
for < i,j < 7. 
In order to prove 3), we must have 



j(O) 



In ® K$(xo,- ■■ ,xr) 



1,(0), 



!,(°) r 



(2/0, 



>!//>(tO-i)®J8 In® Kg (xo, ••• ,asr) 



(?/0,--- ,V„(n)-l)®-fe |i -^Mp(i©j')| = |(iffii%8) -Vi6n((iffii%8)ffi(j©i%8))| 



»(0) 



with y, = xs+i and On = O n , n = 1, 2, 4, 8. The reason 
why we consider this rearranged version is that we show in 



i.e., \(i%8) ■ Tpi6n((i © i%8) © (j © j%8))| = 0. As 8n < 
iffi j < 16n- 1, we have i®j= (16n)(l-2~') +8 l m with 
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In ® La{X0,Xi,X2,Xs) 

In ® ~RJ(X 4 ,X 5 ,X 6 ,X 7 ) 

V -Of(yo, ■■■ ,y p ( n )-i) ®h 



In S> Li(xo,Xi,X2,X 3 ) 
Ol(yo,"- ,Vp(n)-l) <&h 

In ® —Ri(xi, xs, x&, x 7 ) 



In ® i?4(a:4, X5, X6,X7) 

In 8>i4(iBo,a;i,a;a,a;3) 

04n 



Ol(?/0,--- ,l/p(n)-l) ® ^4 ^ 
in ® RJ(X4, X 5 , Xq,X 7 ) 

04n 
ira ® Lj(xo, Xi, X2, X3) ) 



(48) 



Z > 1. So 8 divides i® j as 8 divides both (16n)(l — 2 _i ) and 
8'm. So i%8 = j%8 i.e., * © j = ((» © »%8) © (j © j%8)). 
Thus it is enough to prove that |(i%8) • tpie n (i © j)\ = 0. It 
is indeed true as i(>i6n{i © j) is a multiple of 8. 
To prove 4), we have to show that 

|(8i) • Vi6«(8n © 8» © 8j)| = |(i • V„((i © j). 

We have 8n © 8z © 8j = (16n)(l - 2~ l ) + 8 l m for some Z 
with Z > 1 and m e Z 8 . Let 16n = 2 a and a = 4c + d. 
If Z = c, we have i/j 16 „(8n©8iffi8j) = 8 l X2 d ( rn ) an d ^ n (i® 
j) = 8 l ~ 1 X2 d ( m )- One can easily see that the above statement 
holds. 

On the other hand, if Z < c, we have i/jie n (8n © 8i © 8j) = 
(16n)2- l - 1 + 8 l X 8(m) mdip n (i®j) = n.2~ l + 8 l - 1 xs(m). 
In this case too, the statement holds. 
To prove 5), we have to show that 

\(i ® 8n) ■ ip 16n (i ® j)\ = |»-Vl6n(*©i)|, 

i.e., |(8n) • ipi6n(i © j)| = 0. Now for i © j = or greater 
than 8n, (8n) • ^i6n(i © j) = 0. 
To prove 6), we have to show that 

|0'©8n)-V>16n(«'©j)| = I + \i • 1pl6n(i ® j)\, 

i.e., |(8n)- Vi 6n (»©j)| = 1- But (8n) ■ ipi 6n (i ® j) = 8n for 
all (iffij) e {1,2,3,4,5,6,7,8™}. ■ 

B. Adams-Lax-Phillips Construction from Quaternions and 
Geramita-Pullman Construction as special cases 

Adams-Lax-Phillips has also provided another construction 
of square RODs using Quaternions which is explicitly shown 
inffl. 

Assuming that a square ROD of order n = 2 a 



xo, 



c p(n)-l) 



which has p(n) real variables, is given, then a square ROD 
of order 16n with p(n) + 8 real variables Xi for i = 
0,1,". ,p(n)+7 



16n 



C(^o,- 



i^n)-)-?) 

is given by j48] i, as shown at the top of this page where the 
two matrices L4 and R4 are given by 

/ X Xi x 2 x 3 ^ 

-X x x -x 3 x 2 

-x 2 x 3 x a -Xl 

-X 3 -Xl Xi Xq ) 



L4(x ,xi 1 x 2 ,x 3 ) 



\ 
( 



R4(x4,x 5l xe,x 7 ) = 
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x 5 


x e 


X 7 \ 


x 5 


Xi 


x 7 


-xe 


x 6 


-X-l 


£4 


x 5 


X-j 


xe 


-X5 


X4 j 



respectively with y t = x s+i . 

The Geramita-Pullman construction of square ROD is also 

given by induction explicitly in QJ. 

Consider a recursive construction of square ROD of order 
n = 2 a to 16n as follows: 

Op)=O^)(x ,--,a; p(n) _ 1 ) 



which has p(n) real variables is given, then a square ROD 
of ord 

0, 1, • ■ • , p(n) + 7 is given by 



(GP) 

"i6n °f order 16n with p(n) + 8 real variables Xi for i 



Xg(a;o,- ■■ ,x 7 ) ® I n 



Is * 



,(GP) 



,y P {n)-l) 



Kj{x Q ,- 



(VO, ••• ,2/p(n)_l) 
,Z7) ® /„ 



with 2/j = x s+l . 

It can be checked that both Adams-Lax-Phillips construction 
from Quaternions and Geramita-Pullman's construction given 

. . (Q) (GP) 

in 1 1 1 differ from the constructions of 0^g„ and 16n defined 
above only in rearrangement of variables and in signs of some 
of the rows or columns of the design matrix. 

Lemma 7: Let t > 16 and A t and P t be the square ROD of 
order t as given in Lemma [5] (ii) and (iii), and also let 0^g„ 



and 



^ Gp ) 



be the square RODs given in (l48b and in ( Bit 



»«) 



"16n 



and R 



(GP) 
16n 



which are of order 16n. Then ^4 t 
for t = 16n. 

Proof: Similar to that of Lemma [6] and hence omitted. ■ 

Example 5: Square ROD A\q of size [16, 16, 9] by Adams- 
Lax-Phillips construction from Octonion is given by (l52l . 

Square ROD A 16 of size [16, 16,9] by Adams-Lax-Phillips 
construction from Quaternion is given by (|53V Square ROD 
of _Pi6 size [16, 16, 9] by Geramita-Pullman construction is the 
same as Riq. 

Square ROD of P 32 size [32,32, 10] by Geramita-Pullman 
construction is given by 
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Appendix C 

RATE 1/2 SCALED COD OF SIZE [64, 10, 32] 



17 





■ x 


-£l 


-£2 


-£3 


—£4 


-£5 


-£6 


-£7 


-£8" 


-£16 




XI 


£() 


£3 


-£2 


£5 


— X4 


-£7 


£6 


£9 


£17 




x 2 


-£3 


£0 


£l 


£6 


£7 


— X4 


-£5 


£10 


£18 




x 3 


£2 


-Xl 


£0 


£7 


-XQ 


£5 


— 2:4 


£11 


£19 




X4 


-£5 


-X6 


-XT 


£0 


XI 


£2 


£3 


£12 


£20 




x 5 


X'4 


-£7 


£6 


-Xl 


£0 


-£3 


£2 


£13 


£21 




x 6 


£7 


,X'4 


-£5 


-X2 


£3 


to 


—xi 


£14 


£22 




XT 


-£6 


£5 


X4 


-£3 


-£2 


£1 


£0 


£15 


£23 




Xg, 


-£9- 


-£10" 


-Xll- 


-£12" 


-£13" 


-£14- 


-£15 


£0 


£24 




XQ 


£8" 


-£11 


£10- 


-X13 


£12 


£15" 


-£14 


-XI 


£25 




£10 


£11 


£8 


— £9- 


-X14- 


-£15 


£12 


£13 


-£2 


£26 




m- 


-£10 


£9 


£8" 


-£15 


£14" 


-£13 


£12 


-£3 


£27 




£12 


£13 


£14 


£15 


£8 


— £9- 


-£10" 


-£11 


— 0:4 


£28 




£13" 


-£12 


£15" 


-£14 


£9 


£8 


£11" 


-£10 


-£5 


£29 




X14- 


-£15" 


-£12 


£13 


£10" 


-£ll 


£8 


£9 


-£6 


£30 




£15 


£14" 


-£13" 


-£12 


£11 


£10 


— XQ 


£8 


-XT 


£31 




£i6- 


-£17- 


-£18" 


-£19- 


-£20- 


-£21- 


-£22" 


-£23" 


-£24 


£0 




£17 


£16 


£19" 


-X18 


£21" 


-£20- 


-£23 


£22 


£25 


-£1 




X18" 


-£19 


£16 


X17 


£22 


£23" 


-£20" 


-£21 


£26 


-£2 




£19 


£18- 


-£17 


£16 


£23" 


"£22 


£21" 


-£20 


£27 


-£3 




£20" 


-£21- 


-£22- 


-£23 


£16 


£17 


£18 


£19 


£28 


— X4 




£21 


£20" 


-£23 


£22" 


-X17 


£16" 


-£19 


£18 


£29 


-£5 




£22 


£23 


£20- 


-£21- 


-£18 


£19 


£16" 


-£17 


£30 


-£6 




£23" 


-£22 


£21 


£20" 


-£19- 


-£18 


£17 


£16 


£31 


-£7 




£24" 


-£25" 


-£26" 


-£27" 


-£28" 


-£29" 


-£30- 


-£31 


£16 


-£8 




£25 


£24- 


"£27 


£26" 


"£29 


£28 


£31- 


-£30" 


-£17 


-£9 




^26 


£27 


£24- 


-£25" 


-£30" 


-£31 


£28 


£29" 


-£18" 


-£10 




£27" 


"£26 


£25 


£24" 


-£31 


£30- 


"£29 


£28" 


-£19" 


-£ll 




£28 


£29 


£30 


£31 


£24- 


-£25" 


-£26" 


-£27- 


-£20" 


-£12 




£29" 


-£28 


£31" 


"£30 


£25 


£24 


£27" 


-£26" 


-£21" 


-£13 




£30- 


-£31- 


"£28 


£29 


£26" 


"£27 


£24 


£25" 


-£22" 


-£14 


1 


£31 


£30" 


-£29" 


-£28 


£27 


£26" 


-£25 


£24" 


-£23" 


-£15 


7T 


x 


~ X l 


—x 2 


— £3 


£4 


~ x 5 


-xl 


-x 7 


— x s - 


_:E 16 




X l 


''() 


•'':) 


— x 2 


X 5 


-£4 


— Xj 


;r (i 


' ■) 


X 17 




x 2 


— x 3 


'" 


x l 


;r i 


j; 7 


^4 


~ X l 


;/ 10 


x 18 




£3 


x 2 


~ X l 


■''■() 




~ x l 


x l 


£4 


■'11 


x 19 




xl 


~ x l 


— £ 6 


— £ 7 


■ T 6 


x l 


x 2 


' r :{ 


^12 


'"-•0 




x l 


x 4 


— Xy 


xl 


~ x t 


'■■ 


—x 3 


x 2 


1 x j 


x 21 




x t 


X T 


''4 


~ X l 


—x 2 


•'•':', 


'■''{) 


~ X t 


£ 14 


■' I2 




x* 7 


~ x 6 


X 5 


£4 


~ x 3 


— x 2 


:r\ 


£ 


' T 15 


■" 2 




Xg 


— £g- 


" x 10" 


-£ n - 


~ x \2~ 


_a; 13" 


-£44- 


" X 15 


• r l) 


x 24 




£g 


X S~ 


-£ 11 


x 10~ 


~ x \3 


;/ l^ 


x 15" 


-£14 


~ X l 


I 25 




x 10 


• r ll 


' r r< 


— £9" 


"£44- 


~ x 15 


J; 12 


''-13 


—x 2 


•'^(. 




x 11 - 


" x 10 


• r !) 


£ g - 


" x 15 


x 14" 


" x 13 


''-12 


~ x 3 


x 27 




x \2 


X 'l3 


;/ 14 


X 'l5 


■ T ji 


— Xg~ 


" x 10" 


" X ll 


-£4 


•' ^ 




x \3~ 


~ x 12 


^is" 


-£14 


£9 


• r 8 


■I'll 


_:r 10 


-*8 


x ^.') 




£14- 


_3; 15" 


_a; 12 


x n 


x 10~ 


-£ xl 


;r s 


■ T <) 


-£? 


x 30 




x 15 


x 14 


_a; 13" 


~ x \2 


x \\ 


;/ 10 


~ x 9 


£ g 


-x^ 


''-31 




x 16" 


-x 17 - 


_a; 18" 


~ x l$~ 


~ x 20~ 


" x 21" 


_3; 22- 


" x 23" 


-x 24 


'' r 




x 17 


x 'lfi 


^19" 


~ X 1S 


x 2\~ 


_:E 20" 


_:c 23 


■'■I 2 


' T 25 


" X l 




x 18" 


_:r 19 


;/ 16 


x n 


■ L ±2 


:E 23" 


" x 20" 


~ x 2\ 


'' 26 


— x 2 




x 19 


x i8~ 


-£ 17 


x 'l6 


x 23~ 


_:E 22 


x 21- 


~ x 2§ 


;/ I- 


—x s 




^20" 


~ x 2\~ 


_a; 22- 


~ x 23 


x 16 


;/ 17 


•'-18 


x 'l9 


£28 


£4 




x 21 


x 20~ 


_:E 23 


x 22~ 


-£ 17 


^le" 


_:E 19 


x 18 


; 'i. : ) 


"1 




x 22 


• ; - i -; 


x 20" 


-£21" 


~ X H 


'■' 19 


x 16" 


-x 17 


;/ j° 


x 6 




x 23~ 


~ x 22 


;/ il 


x 20" 


~ x 19~ 


_a; 18 


•'-17 


x 16 


:/ n 


-x 7 




£ 24 - 


_3; 25" 


~ x 26~ 


-£ 27 - 


~ X 2S~ 


_a; 29" 


-£30" 


_x, 31 


;/ i(i 


— Xg 




x 25 


x 24 - 


-x 27 


x 26~ 


-x 2g 


■' if. 


x 31" 


_x, 30" 


-£ 17 


— Xq 




x 26 


x iT 


x 24~ 


" x 25" 


~ X 30~ 


_a; 31 


•' 2'; 


x 29" 


_:E 18" 


" x 10 




x 27 - 


~ x 26 


'■' 2': 


x 24 - 


~ x 31 


^30" 


-x 2g 


x 28" 


_a; 19" 


-x n 




x 28 


•'-!■) 


2 ">M 


''-31 


x 24 - 


_a; 25" 


~ x 2&~ 


-£ 27 - 


_a; 20" 


_:r 12 




x 29" 


~ x 28 


X Z\~ 


~ x 30 


■' l r > 


; 'i4 


x 27 - 


_x, 26" 


" x 21" 


_:r 13 




£30" 


_a; 31~ 


~ x 23 


•'-I-) 


x 2&~ 


-x 27 


J -i4 


x 25" 


_a; 22- 


-£14 




- x 31 


x 30~ 


-x 2g - 


~ x 28 


• r 27 


x 26" 


_:E 25 


x 24" 


_a; 23" 


_:r 15 



